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1 INTRODUCTION 



Laser guide stars (LGSs) l |Foy fc Labeyrie 1985} solve some 
of the drawbacks of conventional adaptive optics (AO) sys- 
tems employing natural guide stars (NGSs), most notably 
sky coverage limitations. However, unlike NGS AO systems, 
LGSs suffer from focal anisoplanatism (FA) or the "cone 
effect" JParenti fc Sasieia 1994ll . This effect becomes more 
severe with increasing telescope aperture size, meaning that 
a single LGS on an extremely large telescope (ELT) would be 
unusable. FA can be at least partially overcome using multi- 
conjugate adaptive optics (MCAO) lIBeckers 198811 . in which 
measurements from multiple LGSs are combined. In this pa- 
per we outline an alternative to a conventional LGS in which 
FA is significantly reduced for a single-conjugate AO sys- 
tem. The technique, which we call SPLASH (Sky Projected 
Laser Array Shack- Hartmann), is a pseudo-reverse of a con- 
ventional Shack-Hartmann wavefront sensor (WFS). Unlike 
a conventional LGS system, the atmosphere is sensed on the 
upward path of the laser. An array of converging beams is 
launched from the primary mirror of the telescope to pro- 
duce an array of Shack-Hartmann spots projected onto the 
sky, which are then imaged by the whole telescope. 

A number of alternatives to "conventional" LGS wave- 
front sensing have been suggested. Baharav et al. 119941 
I199(i|l proposed the creation of a fringe pattern in the at- 
mosphere which is analysed by a Shack-Hartmann wave- 
front sensor. A disadvantage is that very high laser powers 
are required, although this has recently been ameliorated 
by the suggested adaptation of a pyramid wavefront sensor 
l |Ribak fc Ragazzoni 200T| . A proposal by Lloyd-Hart et al. 
involves producing a number of images of different 
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ABSTRACT 

A method for producing a laser guide star wavefront sensor for adaptive optics with 
reduced focal anisoplanatism is presented. A theoretical analysis and numerical sim- 
ulations have been carried out and the results are presented. The technique, named 
SPLASH (Sky-Projected Laser Array Shack-Hartmann), is shown to suffer consider- 
ably less from focal anisoplanatism than a conventional laser guide star system. The 
method is potentially suitable for large telescope apertures ('^ 8 m), and possibly for 
extremely large telescopes. 
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planes in the atmosphere as the laser propagates through 
a focus. These images are then used in a phase diversity 
wavefront sensor. Angel ll2UUlt also proposed a system of 
dynamically refocusing the laser spot so as to effectively in- 
crease the integration time. Kellner et al. (20041 proposed 
using Bessel beams as pseudo-inverse guide stars (PIGS). 
Ribak and Ragazzoni II2004II also proposed a method of re- 
ducing laser spot elongation using distributed launch optics. 
All of these systems share the common characteristic with 
the conventional LGS system that the aberrations are sensed 
during the return downwards path of the laser, unlike the 
scheme proposed here, and they also suffer from focal aniso- 
planatism. 

A further alternative, P** (Projected Pupil Plane 
Patterns) IBuscher et al. 200211 shares similarities with 
SPLASH. A laser beam is expanded to fill the pupil of the 
telescope and is propagated upward through the atmosphere 
as a parallel beam. The beam cross-section is imaged at dif- 
ferent altitudes and the wavefront distortion is determined 
by comparing the intensity distributions. This method does 
not suffer from focal anisoplanatism and, like SPLASH, the 
aberrations are sensed on the upward path. However, 
potentially suffers from two key drawbacks: the observation 
layers need to be separated sufficiently in altitude for the 
beam intensity distribution to evolve, and the measurements 
may be distorted by non-uniform scattering from the atmo- 
sphere. These limitations are not applicable to SPLASH. 

The basic concept of SPLASH was first presented in 
Love et al. 121)041 . but has much in common with Tsch- 
erning aberrometry - a technique used in ophthalmology 
to measure optical aberrations in the human eye. This 
was first described by Tscherning 1189411 and has become 
more widely adopted in recent years, eg. Mrochen et al. 
IMrochen et al. 200011 . Here we present results of a theo- 
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retical analysis and a closed-loop numerical simulation to 
validate the SPLASH technique. In Sect. |5|we outline the 
concept and describe the advantages and problems associ- 
ated with the technique. We present our theoretical analysis 
and its results in Sect.|3| and a description of our numerical 
model and its results in Sect.|l] In Sect.l^we summarise our 
results and outline the remaining problems not addressed in 
this paper. 



There is an additional error on the upward path due to 
a "lever arm effect" . The distance a spot moves due to a lo- 
cal wavefront tilt depends on the altitude at which the tilt is 
applied - a ground layer tilt will cause a larger spot motion 
than a tilt of the same magnitude applied higher up. This is 
not expected to be a significant problem for closed loop op- 
eration, especially for a sodium LGS system, although it will 
have the effect of slightly reducing the system bandwidth. 



2 DESCRIPTION OF SPLASH 

This section describes the key principles of the SPLASH 
technique of wavefront sensing and outlines the advantages 
and disadvantages associated with the technique. 

2.1 SPLASH concept 

SPLASH wavefront sensing is a pseudo-inverse of conven- 
tional LGS/Shack-Hartmann wavefront sensing in which the 
atmosphere is sensed on the upward path of the laser beam 
instead of the return path. The basic premise, illustrated in 
Fig-Q is to project an array of converging laser beams, each 
of size ~ ro (where ro is the Fried parameter) , from the pri- 
mary mirror of the telescope to form an array of spots on 
the sky. The position of each spot on the sky depends on 
the local (subaperture) wavefront gradient. The spots are 
imaged through the full telescope aperture, so the position 
of the final image of each spot will be altered by the global 
(full aperture) tilt. Hence the position of each spot image 
gives a measure of the local tilt minus the global tilt - ex- 
actly the same quantity as is measured in a conventional 
Shack-Hartmann wavefront sensor (WFS) when used with a 
laser beacon. This assumes that the angular size of the spot 
pattern on the sky is smaller than the isokinetic patch (the 
angle over which the wavefront tilt is isoplanatic) . 

The technique could be implemented with either a 
Rayleigh or a sodium beacon, although implementation of 
a sodium layer SPLASH system presents more of a chal- 
lenge in focusing the beams than the Rayleigh version, due 
to the extremely long focal length. Range gating is essen- 
tial, for both the Rayleigh and sodium versions of SPLASH, 
in order to allow the beacons to be imaged without being 
swamped by backscattered light from lower altitudes. 

2.2 Focal anisoplanatism 

The main proposed advantage of a SPLASH LGS is that it 
suffers considerably less from FA than a conventional LGS. 
However, it is still affected by FA and, as a result of the 
unusual configuration, the FA effects are manifested in a 
different way to that seen in a conventional LGS/Shack- 
Hartmann WFS combination. 

SPLASH is affected by FA twice, on the upward and 
downward paths of the light. On the upward path, where 
the wavefront is sensed, the local gradient measured by each 
beam will have a slight error due to FA - see Fig.|21 As each 
spot is then imaged from a different position in the sky, the 
return paths of the light from the spots each sample the 
atmosphere differently so that the global tilt on each sub- 
aperture centroid is separately corrupted by a combination 
of focal and angular anisoplanatism. 



2.3 Effect of diffraction on SPLASH 

Figures and |5| assume geometrical optics but in reality 
the SPLASH beams will be affected by diffraction, result- 
ing in a finite spot size on the sky. The minimum size of 
the spots depends on the size of the launch subapertures, 
the laser wavelength, the laser beam intensity distribution 
and the beacon altitude, with the seeing resulting in a fur- 
ther increase in the size of the spots. Thus the spots could 
potentially overlap on the sky making centroiding difficult 
or impossible. For the sodium beacon case the longer beam 
propagation distance will result in larger sky spots than in 
the Rayleigh beacon case, worsening the problem of overlap- 
ping spots. 

For beams with a Gaussian intensity profile, minimum 
spot size on sky would be achieved by launching the beams 
from subapertures of size ~ 2-3 ro. This would limit the 
wavefront sampling, meaning that SPLASH would be un- 
suitable for high-order wavefront sensing. However, such a 
system would be ideal for adaptive optics in the infrared. 
Launching Gaussian beams would not be possible using the 
approach illustrated in figure^but could be achieved using, 
for example, holographic techniques. 

The problem of overlapping spots could also be reduced 
by implementing a method of time-interleaving the sub- 
apertures such that only a subset of the beams would be 
launched at any given time. All of the spots would be pro- 
jected and imaged in each complete cycle of the system, but 
the wavefront-sensing would be divided into sub-cycles with 
each sub-cycle involving a combination of subapertures cho- 
sen to avoid overlapping of the spots. 

Aside from the issue of overlapping spots, diffraction 
effects may serve to partially remove the effects of FA on 
the upward paths of the beams. The broadened beams will 
sample higher-altitude turbulent layers more fully than the 
idealised geometrical-optics "cones" shown in figures Q and 

m 

2.4 Effect of turbulence on return path 

Return path tip/tilt has already been discussed but some 
consideration needs to be given to higher order aberrations. 
The subapertures in a conventional Shack-Hartmann WFS 
are of similar size to the SPLASH launch subapertures, 
i.e. ~ ro. Thus the wavefront aberration across each sub- 
aperture is dominated by low spatial frequencies and the 
subaperture images are generally not speckled. In SPLASH 
the sky spots are imaged through the full telescope aper- 
ture (with a many-ro diameter) and the turbulence-induced 
aberrations on the downward path can contain significant 
high spatial frequency aberrations. As a result the point 
spread functions (PSFs) could become speckled, and as 
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Spot pattern is reimaged 
at altitude H 



Lenslet array creates spot pattern at 
point optically conjugate to altitude, H 




LASER 
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Returning 
Science Field 



Telescope Primary 



Figure 1. Concept of SPLASH, showing the upward passage of the beams. A possible optical implementation is shown whereby the 
laser is launched via a lenslet array. This is only a conceptual diagram and not a formal optical design - practical implementation may 
require a rather different approach for reasons discussed in the main text. The size of each of the converging beams is ~ ro although we 
have only shown four here for clarity. Furthermore, we have shown the beams as converging to a point, whereas in reality they would be 
broadened by diffraction. See text for more discussion. 



the amount of speckling increases the centroid of the PSF 
can be expected to become poorly correlated with the lo- 
cal wavefront tilt. This effect is known as "centroid aniso- 
planatism" IIYura fc Tavis "19851 IChurnside et al. 19851 . and 
becomes significant for values of D/ro greater than about 10. 

The speckling effect could be reduced by masking off a 
large portion of the aperture so that the spot pattern would 
be imaged through a much smaller aperture. This is not 
really a viable solution, however, because it would lead to an 
increase in tip/tilt anisoplanatism (i.e. differences in global 
tip/tilt on the return paths from the sky spots). 

In a closed-loop system the return-path aberrations 
would be removed by the wavefront corrector, so only the 
residual wavefront error would contribute to centroid aniso- 
plantism. However, if the spots were too badly speckled to 
begin with it might be impossible to close the loop. 



single beacon is split between many WFS subapertures. In 
a SPLASH system, although a separate beacon is launched 
from each subaperture, each beacon is imaged through the 
full telescope aperture and the brightness of the beacons 
can be reduced accordingly. Therefore, providing the beams 
launched from the different subapertures all have the same 
intensity, the increase in collecting area balances the increase 
in the number of beacons. Launch methods in which the 
beam brightness varies between subapertures, such as that 
illustrated in Fig.Q will require higher laser power. 

We now present a theoretical analysis of SPLASH, tak- 
ing into account the effect of FA on the upward path but 
neglecting the effects of diffraction and return path turbu- 
lence. 



2.5 Laser power requirements 

We do not present any formal laser power calculations here, 
however the power requirements of a SPLASH system can be 
expected to be comparable to those of a conventional single- 
LGS system. In a traditional LGS system the light from a 



3 THEORETICAL ESTIMATE OF THE 

EFFECT OF FOCAL ANISOPLANATISM ON 
SPLASH 

This section describes a modal analysis of a SPLASH WFS, 
giving an estimate of the effect of FA on the upward 
(wavefront-sensing) path of the light. The effects of turbu- 
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Figure 2. Left: The upward paths followed by the beams. Each beam samples the atmosphere above its own subaperture, and each beam 
is afEected separately by focal anisoplanatism. Right: The downward paths taken by the light from each spot. The spots are observed 
through the full telescope aperture so the light does not pass through the same section of atmosphere as the upward-propagating 
(wavefront sensing) beams. The upward-propagating beams are affected by the local wavefront gradient, and the downward beams by 
the global gradient (corrupted by focal and angular anisoplanatism). The paths from one subaperture are darkened to show the regions 
of the atmosphere encountered by light from that subaperture. The horizontal lines indicate turbulent layers. 



Figure 3. Illustration of how the telescope pupil is projected onto a turbulent layer at height h = H/A. Areas shaded grey are not 
sensed. Left: Natural guide star - light is parallel so the pupil is fully sampled at all altitudes; Middle: Conventional LGS - as a result 
of FA, the entire pupil is projected onto a smaller circle as altitude increases; Right: SPLASH - each subaperture is projected onto a 
smaller square with increasing altitude as a result of FA, but the spacing of the subapertures remains the same. 



lence on the return path are excluded (but are considered 
later in section |1J. 

Zernike polynomials are a convenient basis set for a 
modal analysis of wavefront correction on a circular aper- 
ture. We will use the same conventions for normalisation 
and numbering of Zernike modes as described by Noll l|197(i|l . 
The phase distortion, ^(JJr), across a circular aperture can 
be expressed in terms of Zernike polynomials, Zj, as 



0(i?r)=^a,Z,(r), (1) 

where R is the radius of the aperture, r is the spatial co- 
ordinate normalised to unit radius and aj are the Zernike 
coefficients, given by 

a, = J ^{Rr)Z,ir)W{r)dv, (2) 
where W{r) is the pupil function. We explicitly exclude the 
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aperture-averaged phase (piston) from the summation in 
equation If the first N Zernike modes could be perfectly 
corrected, the residual wavefront distortion would then be 



JV+l 



(3) 



Providing the Zernike modes are normalised as described by 
Noll II1976II . the mean square residual phase error across the 
aperture can be written as 



2 /J.2 
f^cf, = {<P 



(4) 



A wavefront sensor gives estimates bj of the first A'' Zernike 
coefficients aj . If these modes were corrected as accurately as 
they could be measured, the residual phase variance would 
be 



N oo 

J2{\a',+b',-2a,bA)+ ^"?) 
i=i j=N+i 
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i=i 



JV 

j=i 



(5) 



The variances of the modal coefficients (a^) are given by the 
leading diagonal of Noll's Zernike covariance matrix. Thus 
to predict the performance of the WFS, we need to know the 
variances of the estimated modal coefficients and the co- 
variances between the estimated modal coefficients and the 
actual modal coefficients (ajbj). These variances can be cal- 
culated for a WFS with an NGS using techniques described 
by Wilson & Jenkins lIT^Ml and Cubalchini fl^T^ . 

The phase gradient averaged over subaperture i of a 
Shack-Hartmann WFS is given by 



A 



gi 



ttD ^ 



VZ,(r)dr 



(6) 



subaperture i 



where D is the telescope aperture diameter and A is the 
wavelength (the factor \/tiD scales the phase tilt from units 
of radians of phase per telescope radius to radians of angle) . 
There are two orthogonal phase gradients for each subaper- 
ture (commonly referred to as tip and tilt). This equation 
can be written in matrix form as 

g=4j®°°^ (7) 

where g and a are the vectors of subaperture gradients and 
modal coefficients respectively and is the matrix of {x 
and y) subaperture-averaged derivatives of a large number 
of Zernike functions. Note that, due to the Zernike normal- 
isation used, the elements of are in units of radians of 
phase per telescope radius. The modal coefficient estimates 
b are usually found using 



A 



g 



(8) 



where is the least squares inverse of a i^-matrix con- 
taining a small number of Zernike modes (only as many as 
the WFS is of sufficiently high order to measure). The modal 
covariance matrix for the WFS is given by 



(9) 



where G is a matrix of tip and tilt covariances between 
subapertures (i.e. Gu = (gigt))- Each element of the matrix 
C is a covariance between two Zernike modes, i.e. Cjk ~ 
{bjbk). For a perfect wavefront sensor the matrix C would 
be equal to Noll's matrix (i.e. {bjbk) = {o-jO-k))- 

Up to this point the analysis described has been for 
a Shack-Hartmann WFS viewing an NGS (equivalent to a 
conventional LGS or SPLASH system with turbulence only 
at the ground i.e. with no FA). We are interested in the effect 
of FA on the performance of a SPLASH WFS, so we now 
introduce this into the analysis. We consider an atmosphere 
consisting of just one turbulent layer, although extending 
the analysis to include multiple layers is straightforward. 
We take the height of the layer above the ground to be h 
and the beacon height to be H. Figure |H] shows how FA 
affects the sampling of an atmospheric layer. We define b'j 
to be the SPLASH estimate of Zernike coefficient aj in the 
presence of FA (note that bj still represents the estimate aj 
in the absence of FA). Thus, with reference to Equation |S| 
the residual phase variance is now given by 



oo JV 



JV 



(ajb'j). 



(10) 



We define g[ to be the phase gradient averaged over the 
projection of subaperture i onto the turbulent layer, given 
by 



, A \ ^ 



V^,(r)dr. 



(11) 



subaperture i with FA 



Then g' is the vector of FA-affected subaperture gradients 
given by 



■kD' 



(12) 



where is the matrix of a large number of Zernike func- 
tions averaged over FA-projected subapertures (see Figure 
13 . Two more modal covariance matrices can now be calcu- 
lated: 



where the elements of these G- and C-matrices are 



C' = (^) &'-^G'{&-^f 
C" = {llR^' &'-^G"{^^'-'Y 



(gWi) 



C'jk = (bjb'k) 
Cjk = {b'jb'k). 



(13) 
(14) 



(15) 
(16) 

(17) 
(18) 



We make the assumption that the performance of the 
wavefront sensor is essentially perfect in the absence of FA, 
i.e. flj = bj. Since we are interested in investigating the 
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effect of FA on SPLASH this is a reasonable assumption, 
providing that the number of reconstructed modes is lim- 
ited sufficiently to avoid significant fitting error. Thus the 
diagonals of the matrices C' and C" are equal to the val- 
ues {ajb'j) and (6^^) respectively, which are the unknowns 
required in Eauation llOl 

The matrices G, G' and G" are constructed by map- 
ping tilt covariances for pairs of apertures with appropriate 
spatial separations into square arrays which use the same 
subaperture geometry along both axes as the subaperture 
axis of the ^-matrices. An efficient method of calculating the 
tilt covariance between two spatially separated spatial aper- 
tures has been described by Assemat 1200411 . The tip/tilt 
covariances between subapertures i and I are given by 



{gi,xgi,x) 



1 

2d? 
1 

2d? 



dxdy 



{x,y)(S)Iii{x,y) 



{x,y)(S)Iii{x,y) 



(19) 



(20) 



where the x and y subscripts indicate the direction of the 
tilts. Iii{x, y) is the intercorrelation of the two aperture func- 
tions, defined by 



Iii{x,y) 



u + X V + y 



d 



d 



dudv 



(21) 



where d is the subaperture size and 11(2;, y) is the aperture 
function, given by 



U{x,y) = 1 |a;| < 1/2 and lyj < 1/2 
— otherwise. 



(22) 



D^p is the phase structure function which, assuming Kol- 
mogorov turbulence with an infinitely large outer scale, is 
given by 

D^{x,y) = 6.88r~''/\x' + yY^'^. (23) 

When one or both of the subapertures are affected by FA, 
one of the following equations is used: 



u + x v + y\^( u V . , , 



iii{x,y) 



n 



u + (3x V + f3y 



13d 



13d 



n 



Pd' (3d 



where (3 is an FA factor, defined by 



dwdu(25) 



(26) 



Thus we have all the information required to predict the 
residual wavefront variance for a given number of corrected 
Zernike modes. We have assumed the cone geometry shown 
in Fig.Q and that the system is capable of perfectly correct- 
ing Zernike modes to the degree that they can be sensed. We 
also assume that tip and tilt across the full telescope aper- 
ture can be perfectly sensed, since in reality these modes 
would be sensed using a NGS rather than the laser beacon. 

The results of the theoretical analysis of SPLASH are 
shown in Fig. |1] The benefits of the better sampling of the 
wavefront provided by SPLASH can clearly be seen in the 
plot. The fractional residual variance in each mode (i.e. the 
proportion of the variance in each mode that cannot be 
sensed and corrected) is roughly constant at each spatial 
scale for SPLASH, whereas the effect of FA on a conventional 



LGS becomes more severe as the spatial frequency of the 
aberrations increases. The results indicate that a SPLASH 
system could be expected to perform significantly better 
than an equivalent system using a conventional LGS with 
a Shack-Hartmann WFS, and that the benefits are greater 
for larger telescope apertures. 

It is important to remember that we have ignored the 
effects of atmospheric aberrations on the downward path of 
the light from the focused spots. These effects are included 
in the numerical simulation described in the next section. 



4 CLOSED LOOP SIMULATION 

A closed-loop semi-geometrical Monte Carlo simulation of 
SPLASH has also been implemented. This includes the ef- 
fects of aberrations on the return path of the light through 
the atmosphere in addition to the upward-path FA consid- 
ered in the theoretical analysis. 

For the purposes of the simulation a uniform intensity 
distribution was assumed for the laser beams, with each 
beam being focused on the sky by a lens. The beacons were 
assumed to be sufficiently bright to ignore the effects of pho- 
ton noise and CCD read noise. 

The simulation assumes the diffraction-free geometry 
illustrated in figure|21for the purposes of identifying the sec- 
tions of atmosphere intersected by the beams, but far-field 
diffraction is included in the PSF calculations. For each sky 
spot, and at each time step, the phase aberration induced 
on the upward path is calculated and the on-sky PSF is 
calculated as a FFT of the complex amplitude across the 
subaperture. Similarly, the phase aberration on the down- 
ward path is projected onto the full telescope aperture and 
the downward-propagation PSF is calculated as the FFT of 
the complex amplitude across the telescope aperture. The 
total PSF (the LGS PSF imaged through the atmosphere) 
is calculated by convolving the upward and downward PSFs. 

The simulated spots are centroided separately - they 
are not "stitched" together into a complete spot pattern. It 
is assumed that if the combination of the size and motion of 
the spots is sufficient to cause cross-contamination, this can 
be compensated for by time-interleaving of the spots. 

The wavefront corrector used is an idealised segmented 
deformable mirror (DM) consisting of square segments each 
capable of correcting piston, tip and tilt. Each mirror seg- 
ment is aligned perfectly with a SPLASH subaperture. The 
DM segment tip and tilt values correspond directly to the 
corresponding spot x and y centroids, and the subaperture 
piston values across the DM are reconstructed from the cen- 
troids using a successive over-relaxation (SOR) algorithm 
ISouthwell 198011 . As in the previous section, perfect (NGS) 
global tip/tilt correction is assumed. 

The atmospheric model used consists of two translating 
Kolmogorov phase screens at different altitudes, with 60 per 
cent of the turbulence at the ground and the remaining 40 
per cent at one quarter of the altitude of the laser beacons. 
A range of different D/ro values were simulated and for each 
one, the number of subapertures was matched to D/ro, i-e. 
for D/ro = 8, an 8 x 8 subaperture WFS (and DM) was 
used. 

A conventional LGS AO system with the same DM, 
atmosphere model, beacon altitude, tip/tilt correction, etc. 
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Figure 4. Theoretical prediction of SPLASH performance (solid line) as compared with an equivalent conventional LGS/Shack-Hartmann 
wavefront sensor system (broken line) . Results show residual wavefront variance as a fraction of uncorrected variance in each radial order 
of Zernikes, for a single atmospheric layer at 1/4 of the beacon altitude. Left: 8x8 array of subapertures; Right: 12 X 12 array of 
subapertures. 



was also simulated to allow an objective comparison with 
SPLASH. 

The simulations were carried out on a Cray XDl super- 
computer with 12 Opteron processors running SuSE linux. 
This machine features field programmable gate arrays (FP- 
GAs) for hardware acceleration of numerically intensive pro- 
cessing, but the current version of the simulation code does 
not take advantage of these. Future versions are expected to 
do so IIBasden et al. 2005ll . 

The control loop was successfully closed for all the 
D/ro values simulated and the SPLASH results compared 
favourably with the conventional LGS results, as can be seen 
in Figure |K| SPLASH performs considerably better than a 
conventional LGS system across the whole range of simu- 
lated values. The highest D/ro values simulated correspond 
to an 8 m class telescope, so the technique could be suitable 
for use on existing telescopes. 

As D/ro increases further, the residual wavefront error 
will worsen and centroid anisoplanatism become more no- 
ticeable. For sufflciently large D/ro the residual wavefront 
error will prevent effective wavefront sensing, although the 
regime in which this occurs will be sensitive to the distribu- 
tion of turbulence in the atmosphere. Further simulations to 
higher D/ro values (and with a range of atmosphere mod- 
els) are required to assess the applicability of the SPLASH 
technique to ELTs. 



5 CONCLUSION 

We have described a new method of LGS wavefront sensing 
in which an array of Shack-Hartmann spots are projected 
onto the sky and then imaged through the telescope. 

We have shown theoretically that, in the absence of any 
return-path wavefront aberrations, and assuming purely ge- 



ometrical optics, such a system can be expected to suffer 
considerably less from FA than an equivalent conventional 
LGS system. 

We have further demonstrated the validity of the tech- 
nique using a semi-geometrical closed-loop simulation with 
a realistic atmosphere model, in which return-path aberra- 
tions were included in addition to upward-path turbulence. 
This simulation demonstrates the improvement in perfor- 
mance over a conventional single-LGS AO system for a range 
of D/ro values up to those approximately consistent with an 
8 m class telescope. 

We anticipate conducting further numerical simulation 
work for larger aperture sizes to investigate the feasibility 
of the SPLASH technique for ELTs. It may also be possible 
to experimentally verify the technique on-sky in the coming 
months. 
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